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ABSTRACT 

The very high energy (VHE) 7-ray source HESS J0632+057 has recently been confirmed to be a 7-ray binary. 
The optical counterpart is the Be star MWC 148, and a compact object of unknown nature orbits it every ~321 
d with a high eccentricity of ~0.8. We monitored HESS J0632+057 with the stereoscopic MAGIC telescopes 
from 2010 October to 201 1 March and detected significant VHE 7-ray emission during 201 1 February, when 
the system exhibited an X-ray outburst. We find no 7-ray signal in the other observation periods when the 
system did not show increased X-ray flux. Thus HESS J0632+057 exhibits 7-ray variability on timescales 
of the order of one to two months possibly linked to the X-ray outburst that takes place about 100 days after 
the periastron passage. Furthermore our measurements provide for the first time the 7-ray spectrum down 
to about 140 GeV and indicate no turnover of the spectrum at low energies. We compare the properties of 
HESS J0632+057 with the similar 7-ray binary LS I +61°303, and discuss on the possible origin of the multi- 
wavelength emission of the source. 

Subject headings: binaries: general — gamma rays: general — stars: individual (MWC 148) — X-rays: 
binaries —X-rays: individual (HESS J0632+057) 
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1. INTRODUCTION 

With the advent of the new generation of Imaging At- 
mospheric Cherenkov Telescopes (IACTs) such as MAGIC, 
HESS and VERITAS, a new source class, the 7-ray bina- 
ries, was established. Only few members of this class are 
known to date. Among these objects LS I +61°303, LS 
5039 and PSR B 1259— 63 are regularly detected at very high 
energy (VHE, E > 100 GeV) 7-rays. All of these three 
systems show variable or even periodic VHE 7-ray emis- 
sion, and are spati ally unresolvable by the current gener- 
ation of IACTs dAharonian et al.l J2005t lAlbert et al.ll2006t 
Aharonian et al. 2006; Albert et al. 2009). 



HESS J0632+057 was discovered as an unidentified point- 
like VHE 7-ray source but was considered to be a 7-ray binary 
candidate b ecause of its spatial co i ncidence with the B e star 
MWC 148 dAharonian et al.ll200H Iffinton et alj|2009h . The 
system was observed by VERITAS in VHE 7-rays from 2006 
to 2009 with sparse sam pling and the mea surements did not 
yield any 7-ray signal (lAcciari et al.ll2009l) . The derived flux 
upper limits above 1 TeV were significantly below the pre- 
vious detections, thus suggesting that HESS J0632+057 was 
variable in VHE 7-rays. Since all variable galactic VHE 7-ray 
sources known to date are associated with binary systems^, 
HESS J0632+057 was a very good binary candidate. The 
here in detail presented VHE 7-ray detection simultaneously 
to the X-ray outburst in 2011 Februar y was very recently 
announced by MAGIC and VERITAS llOngl 120111: iMariottil 
l20TltlJogler et al.11201 ltlMaier et alJl201l!K 

Measurements in soft X-rays with XMM-Newton detected 
an X-ray source ( XMMU J063259.3 +054801) at the posi- 
tion of MWC 148 dHinton et alj|2009l) . The X-ray emission 
is well described by a hard power-law spectrum with en- 
ergy spectral index T = 1.26 ± 0.04, consistent with emis- 
sion of synchrotron radiation from VHE electrons, although 
a multi-temperature spectral model can also reasonably de- 
scribe the data. Furthermore, the X-ray source showed a 
variable flux, without changing the spectral shape. A simi- 
lar be havior is seen, e.g., in the 7-ray binary LS I +61°303 
(e.g. lAnderhubetaLll2009l) . Later X -ray observations with 
SwiftfXRT found the source b ut at a different flu x level and 
with a softer spectral index (Falc one et al.ll2010h . Recently 
published Swift/XRT observations from 2009 to 201 1, display 
outbursts in the X-ray light curve fro m HESS J0632+057 with 
a periodicity of P — 321 ± 5 days (iBo ngiorno et al.ll2011l) . 
These measurements also provided evidence for hardness ra- 
tio changes with orbital phase. The periodic X-ray emission 
is a strong evidence for HESS J0632+057 being a 7-ray bi- 
nary. Chandra high time resolution X-ray measurements dur- 
ing the 201 1 February X-ray outburst have been used to search 
for pulsed X-ray emission, but none was detected. Thus the 
nature of the compact com panion of MWC 148 remains un- 
known dRea & Torresll201 11) . 

The region of HESS J0632+057 has also been observed 
at radio wavelengths. The measurements conducted in 2008 
with the Very Large Array (VLA) and the Giant Metre- 
wave Radio Telescope (GMRT) at 5 and 1.28 GHz, respec- 
tively, exhibited an unresolved radio source within the posi- 
tion uncert ainties of the VHE 7-ray source and the Be star 
MWC 148 (Ski lton et alj|2009l) . A flux increase in the 5 GHz 
band from 0.19 ± 0.04 to 0.41 ± 0.04mJy showed the vari- 

33 The Crab Nebula is variable at GeV energies, but no confirmed TeV 
variability has been measured up to now and thus it is not counted among the 
variable VHE 7-ray sources. 



ability of the source on timescales of at least one month. The 
radio data was well described by a power law spectrum with 
energy spectral index a = 0.6 ± 0.2 using non-simultaneous 
data from 1.28 GHz and 5 Ghz. No extended structures be- 
yond the two arcsecond resolution were detected. During 
the 2011 February X-ray outburst very high resolution Eu- 
ropean Very Long Baseline Interferometry Network (EVN) 
observations revealed a point-like source coincident with the 
Be star MWC 148 within uncertainties, which evolved into 
an extended source with a projected size of about 75 AU 
(assuming a 1.5 kpc distance), 30 days later (Mold on et al.l 
I20TT1) . The peak of the emission was displaced 21 AU be- 
tween runs, which is bigger than the orbit size (semi-major 
axis ~ 2.4 AU). The brightness temperature of the source 
was above 2 x 10 6 K hinting to a non thermal origin of the 
particles producing the radio emission. The morphology, size, 
and displacement on AU scales were similar to those found 
in the other gamma-ra y binaries, supportin g a similar nature 
for HESS J0632+057 dMold6netal.ll2011l) . Further high res- 
olution measurements will be needed to understand possible 
morphological changes in the radio structures along with the 
orbital phase. 

Optical radial velocity measurements were taken on 
MWC 148 to verify if it is a mem ber of a binary syste m and 
determine its orbital parameters ( Arag ona et al.l BOlOb . No 
significant radial velocities were found at that time and sim- 
ulations yielded a lower limit on the possible period of the 
system of P > 100 days compatible with the period found 
later in X-rays. Finally, radial velocity measurements with 
the Liverpool telescope obtained from 2008 to 2011 have 
proven the binary nature of HESS J0632+057/MWC 148. 
Fixing the orbital pe riod to 321 days as obtai ned from the X- 
ray meaurements by Bongio rno et al.l (1201 lb these measure- 
ments provide for the first time the orbital parameters of the 
binary system. The compact object orbits MWC 148 on a 
highly eccentric (e = 0.83 ± 0.08) orbit where the perias- 
tron passage occurs at phase (/> pe r = .967 ± 0.008 using 
T = MJD54857.0 dCasares et al.ll2012h . 

In this paper, we present the VHE 7-ray measurements of 
HESS J0632+057 by MAGIC from 2010 October to 2011 
March. In particular we detect VHE 7-rays only during an X- 
ray outburst in 2011 February and measure for the first time 
the spectrum down to 140 GeV. 

2. OBSERVATIONS 

The observations of HESS J0632+057 were performed 
using the MAGIC telescopes on the Canary island of La 
Palma (28.75°N, 17.86°W, 2225 m a.s.l.), from where 
HESS J0632+057 is observable at zenith angles above 22°. 
The MAGIC stereo system consists of two imaging air 
Cherenkov telescopes, each with a 17 m diameter mirror. 
Each telescope features a pixelized photomultiplier tube cam- 
era with a field of view of about 3.5°. The observations were 
carried out in stereo mode, meaning that only shower images 
which simultaneously trigger both telescopes are recorded. 
The stereoscopic observation mode provides a sensitivity so 
that a 5cr signal above 300 GeV is deteced from a source 
which exhibits 0.8% of the Crab Nebula flux in 50 hours ef- 
fective time. The angular resolution is better than 0.07° above 
several hundred GeV and the energy resolution is 16%. Fur- 
ther details on the design an d performance of th e MAGIC 
stereo system can be found in lAleksic et al.l (1201 ll) . 

We observed HESS J0632+057 between 2010 October and 
2011 March for a total of 10.6 hours. All observations were 
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carried out under moonlight conditions and at zenith angles 
from 22 to 50°. The source was observed for several nights in 
each month and each of these observation sets are seperated 
by about 20 days. This strategy maximizes the possibility to 
detect emission from HESS J0632+057 in case of a long or- 
bital period and with the system being active only during a 
short period of its orbit i.e. one observation cycle. Due to bad 
weather no data were recorded in 2010 November and 2011 
January. 

3. DATA ANALYSIS 

The data analysis was performed with the standard MAGIC 
analysis and reconstruction software (MARS). Events that 
triggered both telescopes were recorded and further pro- 
cessed. The recorded shower images were calibrated, cleaned 
and used to calculate image parameters individually for each 
telescope. The energy of each event was then estimated us- 
ing look-up tables generated by Monte Carlo (MC) simulated 
7-ray events. In another step, further parameters, e.g. the 
height of the shower maximum and the impact parameter from 
each telescope, were calculated. The gamma/hadron classifi- 
cations and reconstructions of the incoming direction of the 
primary particl es were perform ed using the Random Forest 
(RF) method dAlbert et al.ll2008l) . The RF calculates the prob- 
ability for each event to be of hadronic origin and denotes this 
parameter as the hadronness of the event. The signal selection 
uses cuts in the hadronness and in the squared angular dis- 
tance between the shower pointing direction and the source 
position (9 2 ). The energy-dependent cut values were deter- 
mined by optimizing them on a sample of events recorded 
from the Crab Nebula under the same zenith angle range and 
similar epochs to HESS J0632+057 data. For the energy spec- 
trum and flux, the effective detector area was estimated by 
applying the same cuts used on the data sample to a sample 
of MC-simulated 7-rays. Finally, the spectrum was unfolded 
in energy, accounting for t he energy resolutio n and possible 
energy reconstruction bias (Albert et al. 20Q7|l. 

The cuts used for producing the 6> 2 -plot for the detection 
were optimized on a Crab Nebula data sample to yield the 
best sensitivity and have a higher energy threshold compared 
to the cuts used to produce the spectrum. For the light curve 
and integral flux calculations we chose a conservative energy 
threshold of E t h = 200 GeV, while the spectrum shows re- 
constructed signals down to 136 GeV. Note that the systematic 
uncertainties at the lowest energies dominate the total mea- 
surement uncertainties. Using a higher energy threshold, 200 
GeV, guarantees smaller systematic uncertainties for the light 
curve and thus yields a better comparison to measurements 
from other instruments. 

4. RESULTS 

We detect VHE 7-ray emission from the HESS J0632+057 
data set recorded in 2011 February, at an orbital phase sep- 
aration of ^0.3 after periastron, with a significance of 6.1a 
in 5.6 hours (see Figure [T). The VHE 7-ray source is not 
resolved by MAGIC and its extension must be smaller than 
the MAGIC point spread function, whose Gaussian sigma is 
0.09 deg above 200 GeV. We obtain an integral flux of F(E > 
200 GeV) = (8.8±1.7 st at±2.1 sys t) x 10~ 12 cm^s" 1 which 
corresponds to abou t 4% of the Crab Neb ula flux. Previously 
reported detections (Aha ronian et al.ll2007l) measured the flux 
only above 1 TeV but agree well within the statistical uncer- 
tainties of our measurements when their spectrum is extended 
to our lower energy threshold. Thus we conclude that dur- 
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FIG. 1. — The squared angular distance between pointing direction of the 
shower and the source position (# 2 -plot) for the position of HESS J0632+057 
(points) and the simultaneous determined background regions (grey shaded 
histogram) for the entir e 201 1 Fe bruary MAGIC data set. The significance is 
calculated according to Li & Ma ( 1983|). N on is the number of events at the 
source position, N g is the number of background events, Af cx is the number 
of excess events (N e x = N on — N a ff ). 

ing our observations, HESS J0632+057 exhibited similarly 
intense VHE e mission to the previous ly detected active VHE 
7-ray episodes d Aharonian et al.ll2007l) . 

The system was only detected in the 2011 February data 
during the X-ray outburst observed by Swift, No indication 
of significant emission was found in the data from 2010 Oc- 
tober, 2010 December or 201 1 March. We denote these three 
months as the non-detection period (NDP). The integration 
time in the individual months of the NDP is, however, rela- 
tively short compared to the 2011 February and we combine 
the NDP to have the highest possible sensitivity for a base- 
line VHE flux. We obtain a flux upper limit for the NDP of 
F(E > 200 GeV) < 3.7 x 10~ 12 cm^s- 1 at th e 95% confi- 
dence level following the method suggested by Rol ke et alJ 
(120051) . Our flux upper limit excludes a baseline emission 
down to the level of 1.7% of the Crab Nebula flux. No in- 
dividual night during the quiescent 7-ray state shows any in- 
dication of a signal. 

We show in Figure |2] the obtained light curve above 200 
GeV of HESS J0632+057 for the nightly averages. The VHE 
gamma-ray source exhibits variability timescales of about one 
month. Faster variability is possible but to detect it a denser 
sampling of the LC is needed. However, no short timescale 
(days) variability is observed during the period of 7-ray activ- 
ity in 2011 February. Under the assumption that there is no 
short-time variability in the time period when no VHE 7-ray 
emission is detected we conclude that the system shows flux 
variations of at least a factor of two between its quiescent and 
active state and that the active state must last between 20 and 
80 days in the VHE regime. 

A correlation of the VHE 7-ray emission with the 2011 
February X-ray outburst is suggestive but can not be proven 
statistically with our sparsely sampled light curve. More ex- 
tensive observations in VHEs are needed for individual night 
correlation studies. Note that only in the time of high X-ray 
activity the system was detected by MAGIC. Whereas the X- 
ray light curve shows a clear peak shape for the outburst, the 
VHE light curve during the 7-ray activity shows a constant 
flux and no variability. Whether this constant 7-ray flux is an 
artifact of the sparse sampling or a real characteristic of the 
outburst cannot be determined with these limited data. Note 
that a similar peak profile as in the X-ray outburst could be 
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FIG. 2.— The light curve of HESS J0632+057 above 200 GeV is shown 
in blue for the MAGIC observations (upper panel) and the Swift/XKT data 
between 2 and 10 keV (lower panel) during the same orbital cycle as the 
MAGIC data in black. Significant emission in VHE is found only in 2011 
February for the duration of about one month contemporary to the X-ray 
outburst. No variability in VHEs is seen during the active state. All error bars 
show the one-sigma statistical uncertainties. 




E [GeV] 



FIG. 3. — Differential energy spectrum of HESS J0632+057 between 
136 GeV and 4 TeV. The horizontal error bars represent the bin width whereas 
the vertical ones show the one-sigma statistical uncertainty. The spectrum is 
well described by the fitted simple power law with T = —2.6 ± 0.3 s iai i 
0.2 sysl shown as the black line. The dashed (blue) line is the spectrum ob- 
tained by H.E.S.S. between 2004 March and 2006 March, and is in good 
agreement with our measurement. 

present in the VHE light curve. 

We obtained a spectrum from the 2011 February data set 
and it is compatible with a simple power law (see Figure |3) 
with photon spectral index 2.6 ± 0.3 sta t ± 0.2 syst and normal- 
ization (1.2 ± 0.3 stat ± 0.2 syst ) • lCT^TcV- 1 



-^s -1 . The 

measured photon index is in very good agreement with the 
one previously published by H.E.S.S. (r = 2.53 ± 0.26 sta t ± 
0.2 syst ), although their spectrum was o btained at energies 
above 400 GeV (lAharonian et all 120071) . No indication of 
a turnover is found in the newly opened energy range by 
MAGIC. 



5. DISCUSSION 

The VHE gamma-ray data obtained during the periodic X- 
ray outburst of HESS J0632+057 that took place during 201 1 
February suggests that the VHE light curve shows similar out- 
bursts like the X-ray light curve. The detection of VHE 7-ray 
emission only during the X-ray outburst suggests a common 
origin although our data are too sparsely sampled to allow 
night to night correlation studies. The MAGIC detection of 
the source during the peak of the X-ray outburst yielded a sim- 
ilar flux level and e nergy spectrum as tho se obtained by HESS 
four years before (Ahar onTan et al.l 120071) . Similar spectral 
shape and flux levels indicate that the same processes might 
be at work during the 7-ray active states. In case of a periodic 
modulation with a period of about 321 days, such a behavior 
would be expected. 

For the first time we could measure the spectrum of 
HESS J0632+057 between 136 and 400 GeV, and find no evi- 
dence for a spectral break or a deviation from a simple power 
law. This is a common feature foun d in other binaries which 
exhibit outbursts (e.g. LS I +61°303 lAlbert et al.ll2009t) . Thus 
the turnover in the spectrum must lie below the energy thresh- 
old of our observations. All currently known 7-ray binaries 
show their maximum emission in the high MeV to GeV en- 
ergy range. This might be the case for HESS J0632+057 
although it has not yet been detected by FermifLPd. The 
integration time required to detect HESS J0632+057 with 
ferm//LAT will depend strongly on the source's duty cycle 
and the spectral properties at MeV to GeV energies. 

In several recent publications the similarity between 
HESS J0632+057 and LS I +61°303 was st ressed based 
on their similar multiwav elength emission (e.g. iHinton et al.l 
120091 iSkilton et al.ll2009h . We note that there might be some 
differences in the VHE emission. LS I +61°303 shows com- 
plex VHE beha vior such as varia bility on timescales as short 
as one day (e.g. lAnderhub et a l. 2009), and different VH E 7- 
ray flux states (lAcciari et al.ll201 It lAleksic et al.ll2012l) . Re- 
vealing similar behavior in HESS J0632+057, requires higher 
sensitivity and better temporal sampled data than available to 
date and thus the degree of similarity between LS I +61° 303 
and HESS J0632+057 might be smaller compared to the con- 
temporary view. Interestingly, the orbital phase lag for the 
detected VHE 7-ray emission in HESS J0632+057 is quite 
close to that in LS I +61°303, about 0.3 after periastron pas- 
sage dLFeF^l] [20TTJ; |Caiare£etal]|20T2). Similar processes 
might produce both the X-ray and the VHE 7-ray radiation. 
However, different spatial- and time-scales of the emitter, 
given the wider and more eccentric orbit in HESS J0632+057, 
and different star-emitter-observer geometries (important in 
the leptonic scenario), make any direct comparison diffi- 
cult. A proper characterization of the radiation and magnetic 
fields, and of possible adiabatic losses, is required. Said that, 
however, we can interpret our HESS J0632+057 data in the 
context of the X-ray data, as done for LS I +61°303 (see 
lAnderhub etai]|2009h . 

In the case of LS I +61°303, the X-ray/VHE emission was 
suggested to originate from a homogeneous leptonic emit- 
ter relatively close to the compact object unde r dominant 
IC los ses. This scenario, already explored in IHinton et al.l 
(2009), is supported by the similar fluxes in X-rays and VHE, 
which is hard to explain in the context of an homogeneous 
emitter producing 7-rays (and the X-ray emitting e -pairs) 
via proton-proton collisions. In the case of HESS J0632+057, 
however, the X-ray luminosity is slightly lower than the VHE 
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luminosity. Although in this case the fluxes and spectra still 
allow the leptonic homogeneous scenario with dominant IC 
losses, a hadronic homogeneous emitter cannot be discarded. 
Therefore the question of the hadronic or leptonic nature of 
the emitter in HESS J0632+057 remains unsolved. 

To distinguish between the hadronic and leptonic pictures, 
and the one population hypothesis for the VHE 7-ray and the 
X-ray emission needs a better sampling of the light curve at 
the relevant orbital phases. A spectral index correlation study 
would shed light on the origin of the emission. Such a study 
will require long and frequent individual night observations 
to provide high enough sensitivity for more precise spectral 
information. Fortunately, due to the X-ray periodicity, such 
measurements can be planned well in advance. Complemen- 
tary information from the GeV band will also be valuable, e.g. 
revealing the spectral shape in the currently unmeasured en- 
ergy range below 100 GeV and above 100 MeV since this in- 
formation could reveal the population(s), of particles involved 



in the emission processes. 



ACKNOWLEDGMENTS 

We would like to thank the Instituto de Astrofisica de Ca- 
narias for the excellent working conditions at the Observato- 
rio del Roque de los Muchachos in La Palma. The support 
of the German BMBF and MPG, the Italian INFN, the Swiss 
National Fund SNF, and the Spanish MICINN is gratefully 
acknowledged. This work was also supported by the Marie 
Curie program, by the CPAN CSD2007-00042 and MultiDark 
CSD2009-00064 projects of the Spanish Consolider-Ingenio 
2010 programme, by grant DO02-353 of the Bulgarian NSF, 
by grant 127740 of the Academy of Finland, by the YIP of 
the Helmholtz Gemeinschaft, by the DFG Cluster of Excel- 
lence "Origin and Structure of the Universe", and by the Pol- 
ish MNiSzW Grant N N203 390834. 

Facilities: MAGIC 



REFERENCES 



Acciari, V. A., et al. 2009, ApJ, 698, L94 

— . 2011, ApJ, 738,3 

Aharonian, E, et al. 2006, A&A, 460, 743 

Aharonian, F. A., et al. 2005, A&A, 442, 1 

— . 2007, A&A, 469, LI 

Albert, J., et al. 2006, Science, 312, 1771 

— . 2007, Nucl. Instrum. Meth., A583, 494 

— . 2008, Nucl. Instrum. Meth., A588, 424 

Albert, J., et al. 2008, ApJ, 674, 1037 

— . 2009, ApJ, 693, 303 

Aleksic, J., et al. 2012, Astroparticle Physics, 35,435-448 
Aleksic, J., et al. 2012, ApJ, 746, 80 
Anderhub, H., et al. 2009, ApJ, 706, L27 

Aragona, C, McSwain, M. V., & De Becker, M. 2010, ApJ, 724, 306 
Bongiorno, S. D., Falcone, A. D., Stroh, M., Holder, J., Skilton, J. L., 

Hinton, J. A., Gehrels, N., & Grube, J. 201 1, ApJ, 737, LI 1 
Casares, J., Ribo, M., Ribas, I., Paredes, J. M., Vilardell, E, & Negueruela, I. 

2012, MNRAS, ArXiv e-prints 1201.1726 



Falcone, A. D., Grube, J., Hinton, J., Holder, J., Maier, G, Mukherjee, R., 

Skilton, J., & Stroh, M. 2010, ApJ, 708, L52 
Hinton, J. A., et al. 2009, ApJ, 690, L101 

Jogler, T., Munar-Adrover, P., Ribo, M., & ParedesJ. M. for the MAGIC 

collaboration. 2011, ArXiv e-prints 1110.1682 
Li, J.,etal. 2011, ApJ, 733, 89 
Li, T.-P., & Ma, Y.-Q. 1983, ApJ, 272, 317 

Maier, G, for the VERITAS Collaboration, Skilton, J., & for the HESS 

Collaboration. 201 1, ArXiv e-prints 1111.2155 
Mariotti, M. 2011, The Astronomer's Telegram, 3161, 1 
Moldon, J., Ribo, M., & Paredes, J. M. 2011, A&A, 533, L7 
Ong, R. A. 2011, The Astronomer's Telegram, 3153, 1 
Rea, N, & Torres, D. F. 201 1, ApJ, 737, L12 

Rolke, W., Lopez, A., & Conrad, J. 2005, Nucl. Instrum. Meth., A551, 493 
Skilton, J. L., et al. 2009, MNRAS, 399, 317 



